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Ph ' Abstract. We report very preliminary results of a Dalitz-plot analysis [1] of A+ in the decay to 

\ pK~K + with the helicity formalism. We used the data from the fixed target experiment SELEX [2] 

(E78 1) in Fermilab. We report about branching-ratios of the resonant states involved, and a possible 
— ■ initial state polarization. 
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INTRODUCTION AND ANALYSIS METHOD 

(N : 

^ ! In our search for resonances in the 3-body decay of A+ — > pK~n + we are considering 

in addition to the non-resonant mode the following resonant decay modes: A+ — > K*°p, 
K*° -> K A+ -> A++K , A++ -> pn+; and A+ -> A(1520)ji+ A(1520) -> p^. 

V ' As a search tool for the resonances we calculate the invariant masses of pairs of 

daughter particles M?-, Atf k and fill with them into a two-dimensional histogram (Dalitz 

r-| ! plot). A Monte-Carlo simulation of what could be expected is shown in fig. 1 (right). 

The data sample used for this analysis is the same as in a previous SELEX publica- 
tion [3], where more details can be found. We applied the following cuts to get a clean 
A+ signal: good fits for tracks and vertices; momentum > 8.0GeV/c for all tracks; pro- 
ton and kaon identified in the RICH [4]; secondary vertex outside material; separation L 
of primary and secondary vertices L > 8 o, where o denotes the combined error of the 
vertices; a < 1.7 mm; at least two tracks with a miss distance to the primary vertex of 
more than 20 /jm; the momentum vector of the reconstructed A+ has to point back to 
the primary vertex. The obtained pK~n + invariant mass distribution can be found in 
fig. 1 (left). 

To eliminate the contribution of background under the A+-peak, we produce a Dalitz- 
plot of the sidebands shown in fig. 1, with a proper mapping of the allowed phase space. 
After normalizing to the correct number of events, we subtract this contribution bin-by- 
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FIGURE 1. Left: Invariant mass distribution for pK~K + . A signal of approximately 1500 A+ can be 
clearly seen. The signal and the sideband regions used in this analysis are indicated. Right: Dalitz plot 
from Monte Carlo with a mixture of non-resonant and three resonant decay modes. 



bin from a Dalitz plot of the signal region. We verified this procedure by comparing 
the plots obtained by different sideband regions on the left and right side of the signal 
region. In fig. 2 we show the Dalitz plots for the sideband and the signal regions. In the 
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FIGURE 2. Left: Dalitz plot for the sideband region indicated in fig. 1 with phase space mapping. Right: 
Dalitz plot of signal region after sideband subtraction. 

background subtracted signal region one can still see some enhancement on the lower 
left of the phase space, also present in the sideband; we are still investigating the nature 
of this. 



To extract the information from the Dalitz plot, we used an helicity formalism [5, 6]. 
We will describe the fit functions used in the following sections. 



HELICITY FORMALISM FOR TWO BODY DECAYS 

From Fermi's Golden Rule the partial decay width is given by dT = ^^-\Q.\ 2 d^> n , with 
dT ~ |Q| 2 = | <BC\T\A > | 2 , where A and BC denote the initial and final states. We are 
working in the eigenstate base \M,S,p fJ ,X> of the operators M 2 = P M P M , S = —a> M (if 
with G) a = i£°^ vA - Mjj V P x , and P^CQ = J X P X + J y P y + J Z P Z ; (0° represents the helicity 
operator with eigenvalue X. For simplicity we assume only X>. 

In the rest frame of the mother particle A we can express the spin states as | ja»M > 
and the final state as \Q B § B X B Xc >> with 

E V -^t^mx, (to, e s , -$b) \jmx b x c > . 

JM V 47C 

Applying angular momentum conservation, we can write the transition amplitude as 

<5C|r|A>= ^y±^D*J* ki (to,Q B ,-to) <X B X c \T\j A m A > 

With D^((|),0,-(|)) = e-^ m -Vd J mX (Q), Xi = X B -X C (M < m A ), and summing over 
all spin and helicity projections and we obtain 



HELICITY FORMALISM FOR THREE BODY RESONANT 

DECAYS 

The decay width in this case is given by dT ~ | < DE\T2\B >< BC\T\ \A > | 2 . In the rest 
frame of the resonance B, with the z-axes pointing into the direction of motion of the 
resonance B in the rest system of the mother particle A, we can write 

<DE\T 2 \B>= y^5±i e ^(^-^)^ 2 (0^) <X D X E \T 2 \B> . 
Summing over all resonances we obtain 

m A X(jXj)Xe Xg B 

(1) 

where we also consider the finite width of the resonance with a Breit-Wigner distribution 
BW(m r ) ~ , m f r , — — . In the parity conserving decay of the resonance we can write 

ax D x E = (-l) SD+SE - SB T\Bn D T\Ea_x D _x E . 



In this formalism we can naturally introduce an initial polarization of the mother 
particle P A , given by dT ~ \{\ +P A ) ^ XdXe \ L^BW^^i^^/ 

+2(1 - p A)LXcX D X E \Lx B LBBW(mB)^ By _i jXB ^ XDjX f. 

Integrating the contribution of a resonance over the phase space gives 



F r = 





BW{m r )^ mA x B 


2 dx 


I ^niA ,~Kb 


LBBW(m B )^ B;inA x B \ 2 dx 



(2) 



from where, after applying weight factors for isospin conservation, we can extract the 
branching ratio for the resonance. 



PRELIMINARY RESULTS 



We performed an unbinned maximum likelihood fit with the functions described in the 
previous section, with 23 free parameters in the fit. The preliminary results are shown in 
the following tables. Only statistical errors are shown. 





Parameter 


Amplitude 


Phase 




Ni 


1. fixed 


0. fixed 


nonres. 


N 2 


320 ±82 


.2±1.9 




N 3 


26 ±56 


4.7±1.9 




N 4 


200 ±135 


3.2±0.2 


K*° 


Ax 


495 ±143 


5.6±0.2 




A 2 


70±91 


3.9 ±10 3 


K*° 


A3 


362 ±91 


3.3 ±0.1 




A 4 


95 ±69 


2.9±481 




Bi 


11 ±23 


3.3 ±10 4 


A++ 


B 2 


196 ±87 


3.5 ±.9 




Ci 


115 ±142 


2.6±18 


A(1520) 


c 2 


644 ±144 


0.8 ±.2 



Polarization P A = 0. 1 ± 0.4 



Branching Ratio 


nonres. 


0.73 ±0.29 


K*° 


0.14±0.17 


A++ 


0.09 ±0.14 


A(1520) 


0.04 ±0.09 



The two most significant features in the Dalitz plot are not properly taken into account 
in this analysis method. The fit function does not include a resonance describing the 
feature seen in the lower left of the phase space, leading to an over-estimate of the 
non-resonant contribution, and the central region with a small number of entries is not 
taken into account in an unbinned fit. At this moment, we are finalizing the analysis 
procedure, by optimizing the cuts used for extracting the A+ signal, a study if additional 
resonances have to be included, and on a binned fit to included the information about the 
center region. Also some systematic studies are under way to quantify the significance 
of the results. This work was supported by CONACyT Mexico under Grant 28435-E. 
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